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The Josephson effect is theoretically studied in two types of SQUIDs consisting of s wave super- 
conductor and Sr 2 Ru04. Results show various response of the critical Josephson current to applied 
magnetic fields depending on the type of SQUID and on the pairing symmetries. In the case of a 
p x + ipy wave symmetry, the critical current in a corner SQUID becomes an asymmetric function 
of magnetic fields near the critical temperatures. Our results well explain a recent experimental 
finding [Nelson et. al, Science 306, 1151 (2004)]. We also discuss effects of chiral domains on the 
critical current. 
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I. INTRODUCTION 

One of the important developments for unconventional 
superconductivity after the the discovery of high-T c su- 
perconductivity, has been the series of so-called phase 
sensitive experiments which test for the symmetry of the 
Cooper pair wave function. These are the SQUID-type of 
interference experiments^! 4 ^, the observation of sponta- 
neous half-flux quantization in frustrated loops^ and the 
measurement of zero-bias peaks in quasiparticle tunnel- 
ing spectra indicating subgap quasiparticle states at the 
sample surfacei2i2iiSiii. This set of experiments, techni- 
cally rather diverse, is based on the same concept, the un- 
conventional phase structure of the superconducting con- 
densate, and have uniquely proven that the Cooper pairs 
have the spin-singlet d x i_ y i wave symmetry. A further 
unconventional superconductor whose pairing symmetry 
has been established with high confidence is S^RuO^i 
This is a spin triplet p-wave superconductoBiiiiii with 
a gap function of the form d(k) oc z(p x ± ip y )i£iiL^ a 
so-called chiral p-wave state with whose Cooper pairs 
possess an angular momentum component along the z- 
axis. Tunneling experiments have shown the presence 
of subgap surface statesA2*A2i2Sl. Moreover, the anoma- 
lous temperature dependence of the critical current in 
Pb/Sr2Ru04/Pb Josephson junction arrangement 21 has 
been interpreted in terms of an interference effec1i22iS2iS4. 
A direct experiment of the type of a SQUID-interference 
or frustrated loop is difficult here for many technical rea- 
sons and has not been performed until very recentl}i2&. 

In this paper we would like to analyze some issues 
which have to be taken into account in the interpreta- 
tion of the SQUID-type experiments for Sr2Ru04. The 



basic principle had been designed long ago by Geshken- 
bein and co-workers 26 . One of the problems lies in the 
Josephson junctions between a conventional spin-singlet 
s-wave superconductor and a spin-triplet p-wave super- 
conductor. The mismatch of the angular moment (par- 
ity) and of spin quantum number on the two sides of a 
Josephson junction seems at first sight inhibit the lowest 
order Josephson effect so that only a coupling in second 
order would be allowed. This would indeed be fatal for 
phase sensitive tests based on the Josephson effect. It 
has, however, been shown that the presence of spin-orbit 
coupling saves the situation since only the total angular 
momentum has to be conserved. Thus under well-defined 
conditions, the lowest order between an s- and a p-wave 
superconductor is possibl o 27 i 28 i 29 i 30 i 31 . The conditions 
leave a certain arbitrariness concerning the sign of the 
Josephson coupling which can be important for interfer- 
ence effects. This can be illustrated, if we consider the 
definition of the lowest order matrix element which can 
be derived from a simple microscopic tunneling model: 

(Mk)(kxn)-d(k)) FS (1) 

where the average runs over the Fermi surface and n is 
the interface normal vector. The arbitrariness appears in 
the orientation of the normal vector, into or out of the 
p-wave superconductor. This has been recently demon- 
strated by Asano and co-workers 31 for a model where the 
spin-orbit coupling of the interface potential was taken 
into account for this matrix element whose sign then de- 
pends on the shape of the interface potential. Thus de- 
tails of the spin-orbit coupling and the interface poten- 
tial, e.g. the way the parity is broken at the interface, 
matter for the Josephson phase relation. The arrange- 
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ment suggested by Geshkenbein et al. relies on the as- 
sumption that all interfaces treat parity the same way 26 . 
We will follow this assumption here too. 

We discuss now two basic forms of SQUID-interference 
devices built from an s- and a p-wave superconductor. 
The first type in Fig. [IJa) is the symmetric SQUID as 
proposed by Geshkenbein and co-workers 26 and the sec- 
ond in Fig.^b) is the corner-SQUID analogous to the one 
used for high- T c -super conductors. The p x + ip^-phase of 
Sr2Ru04 introduces an additional problem. This phase 
is two-fold degenerate, i.e. p x + ip y and p x — ip y are 
equivalent and would in general form domains in a sam- 
ple, depending on the cooling history. We will discuss in 
the following also the implications domains on the inter- 
ference experiment. Similar problems appear for other 
chiral superconducting phase as we will discuss below. 

This paper is organized as follows. In Sec. II, the 
Josephson effect in two types of SQUID is discussed for 
Px + Wy-> Px and p y pairing symmetries. Effects of the 
chiral domain are studied in Sec. III. In Sec. IV, we dis- 
cuss the critical current in the chiral d and the chiral / 
wave symmetries. The discussion and the conclusion are 
given in Sec. V and XI, respectively. 



II. JOSEPHSON CURRENT 

Before discussing the SQUID-experiment, the basics 
of the Josephson current-phase relation between s wave 
superconductor and S^RuCU (SRO) in Fig. [He) has to 
be addressed. We describe the gap function of Sr2Ru04 

A p =i(rj x p x + rjyp y )z ■ &a 2 (2) 
=iA(p x ± ip y )z • <t<j 2 (3) 
=iAe ie z-aa 2 , (4) 

where aj with j = 1,2 and 3 are the Pauli matrices, rj x 
and rj y are the two complex order parameters and z is 
taken as a unit vector parallel to the z-axis™ Moreover, 
Px = Px/pf — cos 6 (p y = Py/pF = sin#) is the normal- 
ized momentum component on the Fermi surface in the 
x (y) direction with pp being the Fermi momentum of 
the p-wave superconductor. Assuming a cylindrical sym- 
metric Fermi surface, we can represent this gap function 
also simply to the angle on the Fermi surface, reflecting 
best its internal phase structure. The gap function of the 
5-wave superconductor is given by = iA<72. Without 
any loss of generality, we may take the gap magnitudes 
identical in both superconductors, A. On the basis of 
the current-phase relationS2i22i^, the Josephson current 
for the lowest two orders derived from a microscopic cal- 



culation close to T c 31 can be written as 

Jp x ±ip y =Ji cos (^ + On) ~ J2 sin (2((f + n )) , (5) 

jl = TB as (^) t> (6) 
J2=T *(^) 3 ^' (7) 

where 6 n is the angle of the junction normal vector in the 
plane relative to the x-axis, Tb denotes the transmission 
probability of a Cooper pair and as is a measure for 
the strength of the spin orbit coupling. Note that it can 
have either sign depending on the junction. The Joseph- 
son current J is measured in units of eAoN c /h : where 
Ao is the amplitude of the gap function at T = and 
N c is the number of propagating channels on the Fermi 
surface. Constant coefficients of the order of unity have 
been omitted in Eq. (0. From Eq. it is clear that 
the coupling between s-wave and p x + ip y -wdive super- 
conductor for a junction in x-direction (0 n = 0) goes via 
the p y component. We assume a gauge where the overall 
phase of the p-wave order parameter is shifted by ±tt/2 
for the p^-component. Thus, in lowest order the current- 
phase relation is proportional to coscp rather than sin(p, 
in Eq. (5J). 

In the first step, we consider the interference pattern 
of the two basic SQUID's in Fig. [Ifa) and (b). We as- 
sume that all junctions are of the same in the sense, that 
the normal vector n in Eq. is pointing a definite di- 
rection, say towards the s-wave superconductor. In the 
symmetric SQUID, the Josephson current of the left and 
right junction are given by 

Jl(^) =Ji cos ip — J2 sin 2ip, (8) 
Jr(v) — — Ji cos cp — J 2 sin 2(p. (9) 

The relative sign change between the first-order terms on 
both sides is due to the the corresponding n which differ 
by 7r. The Josephson current in the symmetric SQUID is 
then expressed 

J s (<p, *) + (^b) + Jr(<P ~ ^b), (10) 

= — 2 Ji sin cp sin (\>b — 2 J2 sin 2cp cos 2(/>#, (11) 

where <pB — tt&/&o an d $0 = 2ithc/e. In the corner 
SQUID, the Josephson current in the bottom junctions 
is given by 

Jb(p) = — Ji sin<£ + J2 sin2<p, (12) 

where the current-phase relation is derived from Eq. @ 
with n = 7r/2. The Josephson current in the corner 
SQUID is expressed in the same way, 

Jc(<P, =Ji cos((^ + 0b) - Ji sm(cp - 0b) 

— 2 J2 cos 2(/?sin20 j B. (13) 

Note that for the corner SQUID it matters which state is 
realized. For Eq. ((12)) we assumed that the state p x +ip y . 
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in high temperatures near T c . Since J\ ^> J2, we find, 



(a) Symmetric SQUID 
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(b) Corner SQUID 
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FIG. 1: Schematic pictures of two types Josephson junction 
are shown in (a) and (b). The singly connected Josephson 
junction is given in (c). In (d)-(f), the SQUID of p x and p y 
symmetries are shown. The Josephson effect in the corner 
SQUID of p y symmetry is identical to that of the p x symme- 
try. 
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FIG. 2: The Josephson critical current is plotted as a function 
of $ in the symmetric SQUID in (a) with J1/J2 = 10. The 
calculated results for the corner SQUID are shown in (b). 



while for p x — ip y would yield a sign change of the first 
term. 

In Fig. [21 the critical Josephson current in the sym- 
metric SQUID and that in the corner SQUID are shown 
as a function of A parameter J1/J2 = 10 is realized 



max|J s ($)| ^2Ji 
max|J c ($)| ^2Ji 



sin 7T 



$0 



sin ( 7r— - ± — 

$0 4 



(14) 
(15) 



The odd parity symmetry immediately results in a min- 
imum of the critical current at $ = in the symmet- 
ric SQUID2&, which is connected with the opposite sign 
of Josephson coupling on the two junctions according to 
Eq. (PJ. Note that this pattern is not dependent on which 
of the two degenerate p-wave state is realized. Actually 
for the symmetric SQUID any p-wave pairing state gives 
the same qualitative behavior and does not depend on the 
detailed symmetry, such as broken time reversal symme- 
try, as long as the first-order Josephson coupling induced 
by spin-orbit effects is dominant. This is different for 
the corner SQUID which leads to a distinction between 
different p-wave states. 

The change to the corner SQUID configuration yields a 
phase shift by ±7r/2 for p x ± ip y . As a consequence, the 
coscp current-phase relation can be exchanged by sincp 
for the bottom junction. For comparison, the results for 
the p x - and p^-wave states are depicted in Fig.[2]with the 
corresponding device illustrations in Fig. E{d)-(f). The 
Josephson current-phase relation in the single junctions 
can for fixed gauge be given by2L24 



4 (<p) 



- — J2 sin 2<p, 

-Ji sin if — J2 sin 2ip. 



(16) 
(17) 



For the p^-component (n parallel to x), there is no con- 
tribution due to spin-orbit coupling, in contrast to the p y - 
component which has such a term proportional to sin (p. 
in the p y symmetry. 

In the symmetric SQUID of the p^-type, the critical 
current has the period of $o/2 as shown in Fig. EJa) 
because there is only the second-order contribution in 
the Josephson coupling. On the other hand, the p y -type 
behaves identical to the p x + ip y -type. For the corner 
SQUID configuration the ^-symmetry (and equivalently 
for the ^-symmetry) has again a <£>o/2-periodic interfer- 
ence pattern of the critical current. 

If we assume that J\ <C J2 due to weak spin-orbit cou- 
pling, the features of the interference pattern are signifi- 
cantly modified. In Fig. |3 the critical Josephson current 
is plotted as a function of $ for J\ = O.IJ2. For J\ <C J2, 
we actually find 



max|J s ($)| ~2J 2 
max|J c ($)| ~2J 2 



cos ( 2tt — 
$0 

sin ( 2tt — 
$0 



(18) 
(19) 



for the p x + ip y symmetry. The period of J c is given 
by $0/2 because the Josephson current proportional to 
sin 2ip is dominant. The critical current in the symmetric 
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FIG. 3: The Josephson critical current is plotted as a function 
of $ for Ji/J 2 = 0.1. 



SQUID takes its maxima at $ = because sin 2ip remains 
unchanged under the tt phase shift. On the other hand in 
the corner SQUID, the critical current takes its minima 
at $ = 0. This is due to the sign change of sm2cp under 
the tt/2 phase shift in the p x J rip y symmetry. The critical 
current in the p x and p y symmetries can be described by 
Eq. (THJ) irrespective of types of SQUID as shown in the 
Fig. 03 Thus the minima at $ = in the corner SQUID 
directly suggests the p x -r-ip y symmetry in SRO when the 
period of the oscillations is $>o/2. 
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domains on the Josephson current. We consider a simple 
model of such chiral domains as shown in Fig.[SJ where f 
and I indicate domains of the different chiral states. The 
size of domains should be much larger than the thickness 
of domain walls which is typically given by the coher- 
ence length of superconductors. The structure of domain 
walls has been investigated based on Ginzburg-Landau 
theories^. For a domain wall with normal vector n, the 
p-wave component parallel to n, p\\ , keeps its phase, while 
the component perpendicular, p±, changes the phase by 
tt. This has important implication to the analysis of the 
above SQUID configurations which we will discuss here 
for a few simple situations as shown in Fig. [51 under the 
assumption J\ ^> J2. 
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FIG. 4: The Josephson critical current is plotted as a function 
of $ for p x + ip y . 

In Fig.QJ the critical current for the p x -r-ip y symmetry 
is shown for several choices of J1/J2. In the vertical axis, 
successive plots have been offset. The critical current is 
always symmetric with respect to in the symmetric 
SQUID. The current maxima at $0 for J1/J2 = 0.1 is 
changed to the minima as increase of J\ / J2 . In the corner 
SQUID, the asymmetry in the critical current gradually 
disappears with decreasing of J\ j J2 . 

III. CHIRAL DOMAINS 

We have assumed so far that SRO is a single domain 
of the p x + ip y symmetry. Real materials, however, may 
have multi-domain structures of the p x + ip y and the 
Px — Wy symmetries. Here we discuss effects of chiral 



FIG. 5: A simple model of chiral domains is illustrated in the 
symmetric SQUID (a) and the corner SQUID (b). Across the 
domain wall with a normal vector n, the p-wave component 
parallel to n keeps its phase, while the component perpendic- 
ular to n changes the phase by tt. 



The presence of horizontal domain walls in Fig. i.e. 
n || y , implies that the two chiral domains have either 
Px + i>Py or — p x + ip y form. According to our previous 
discussion, no qualitative change of the interference pat- 
tern is expected. Thus the critical current in Fig. 2(a) 
remains unchanged even in the presence of the horizon- 
tal domain walls in the symmetric SQUID. On the other 
hand, the vertical domain walls give rise to domains with 
Px + ipy and p x — ip y form. In this case the number for 
domain walls between the two junctions matters since 
each domain wall switches the phase of p yi entering the 
first-order term of the Josephson coupling, by tt. Thus, 
the interference pattern would look different for an odd 
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FIG. 6: The critical current in the presence of domain walls 
for J1/J2 — 10. In the symmetric SQUID, the results for 
the chiral p wave symmetry in the presence of the vertical 
domain walls are shown in (a). In (b), the critical current for 
the chiral p wave symmetry in the corner SQUID is shown. 
The solid (broken) line denote the results when the chirality 
of the last domain at the bottom is p x + ip y (—p x + ip y ) in 
(b). 



or even number of domain walls: 



max|J s ($)| 



2Ji 
2Ji 



Ndw even 



sin (tt^) 
cos (tt^) N dw odd, 



(20) 



where Ndw is the number of the vertical domain walls. If 
the number of walls is an even integer, the critical current 
takes its minimum at $ = as shown in Fig.EJa). How- 
ever, the critical current has its maximum at $ = when 
the number of domain walls is an odd integer. Thus the 
existence of the vertical domain walls can change drasti- 
cally the characteristic feature of the critical current in 
the symmetric SQUID. 

In the corner SQUID, the current-phase relation in the 
two junctions as shown in Fig. Efb) are given by 

Jl(<p) —J\ cosy? — J 2 sin2(/?, (21) 
Jb(v>) — =t J\ sin<p + J 2 sin2(f. (22) 

where the sign of the first-order term of the bottom junc- 
tion is determined by the chirality of the last domain at 
the bottom. The number of domain walls is irrelevant 
here. Still the asymmetry of the interference pattern in- 
dicates the broken time reversal symmetry in any case as 
shown in Fig. Ob). 



IV. OTHER CHIRAL SUPERCONDUCTING 
PHASES 

The chirality of superconductivity reflects the internal 
angular momenta of a Cooper pair. Also other super- 
conducting phases besides the chiral p-wave phase can 
be found with this property. In this section, we compare 
the Josephson effect of the p x +ipy-wave phase with those 
of other chiral phase of d- and /-wave origin. 

A chiral phase in the case of spin-singlet d-wave sym- 
metry can be composed of the d x i_ y i and the d^-state, 



when they are degenerate, yielding a d x i 
gap function was discussed as the surface states of high- 
T c cuparates 36 where the degeneracy is not given, but one 
of the components is subdominant. The other example 
is Na^CoC^ -^H^O 37 . Thus we consider the gap function 



= {Vx(pl - P 2 y ) + T] y 2p x p y } icr 2 

--A(pl -p 2 y ± i2p x p y )ia 2 = Ae ±i29 i& 2 . 



(23) 
(24) 



Here the coupling to the s-wave superconductor is sim- 
pler, since it is not relying on spin-orbit coupling. Thus, 
the current phase relation dominated by the first-order 
coupling. 



J(4>) = J d sin(0 + 2<9 n ), 
A 



J(i = Tb 



IT) 



(25) 
(26) 



Concerning the domain wall structure of this state, the 
two order parameter components behave in the same way 
as in the p-wave case. Thus only the component perpen- 
dicular to the domain wall normal vector changes sign. 
Thus, for n || x the component rj y (d xy ) changes sign 
while d x 2_ y 2 keeps the phase. Thus, irrespective of the 
domain wall number for the vertical domain wall case, 
the maximum of the critical current in the symmetric 
SQUID would be at $ = 0. Thus, domain walls could 
not lead to an interference pattern like in the p-wave 
case. In the case of horizontal domain walls, the situa- 
tion is more complicated. There is a staggered phase for 
the d 



x* — y* 



component. This would give rise to compen- 
sating contributions of the two domains to the Josephson 
effect, similar to the situation discussed in the context of 
c-axis junction between an s-wave and a twinned high- 
temperature superconductor^^. In addition, this con- 
figuration could give rise to spontaneous half-quanta flux 
lines wherever a domain wall hits the junction perpendic- 
ularly. Under these circumstances, the interpretation of 
interference pattern would require much more care. Note 
that such spontaneous fluxes do not occur in the p-wave 
case as long as the domain wall ends perpendicularly on 
the junction. In both cases, the flux magnitude depends 
on the angle between junction interface and domain wall. 
A similar problem with the domain walls occurs for the 
corner junction as a simple examination reveals. 

The spin-triplet chiral /-wave symmetry is also a possi- 
ble candidate for the superconducting phase in Na^CoC^- 
2/H2O. The proposed gap function is given by 



Af) 



~~ iA {(pI ~ 3p 2 y )p x ± i(p 2 y - 3pl)p y } z • aa 2 
--Ae Tise iz-aa 2 . (27) 



The current phase relation for the chiral / wave symme- 
try is given by Eq. (0. In both the symmetric and the 
corner SQUID, the characteristic behavior of the criti- 
cal current for the chiral p wave symmetry discussed in 
Sec. II and III are also valid for the chiral / wave sym- 
metry. In the limit of J\j J 2 ^> 1, the critical current in 
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the symmetric SQUID take the minima at <I> = and 
that in the corner SQUID shows the phase shift by 7r/4 
in the absence of the domain walls. In addition, there 
are no differences in effects of domain walls on the crit- 
ical current between the chiral p wave and the chiral / 
one. To distinguish the chiral p wave symmetry from the 
chiral / wave, we should design SQUID's with a = 7r/3 
or 27r/3 40 . The critical current for / waves becomes sym- 
metric function of <I> and takes its minimum (maximum) 
at $ = in the SQUID with a = tt/3 (a = 2tt/3) 40 as 
shown in Fig. [7| On the other hand, the critical current 
for J1/J2 ^> 1 in the chiral p wave results in 



max|J a=7r/3 (<I>)| ~ 2Ji cos I tt— - - 
max|j a=27r/3 (<I>)| ~ 2Ji 



$ 7T 
$0 " 3 



(28) 
(29) 



The calculated results are shown in Fig. [7| The asym- 
metry in the critical current persists in the chiral p wave 
in such SQUID's. This argument, however, is valid when 
the difference in the transmission probabilities of the two 
tunnel junctions in Fig. [7| are much smaller than them- 
selves. 




FIG. 7: The critical current for the chiral p wave is compared 
with that for the chiral / wave in (a) and (b), where J1/J2 = 
10. The relative junction angle is set at a = ty/3 in (a) and 
a = 2tv/3 in (b) as shown in the upper figures. 



around $ = in the symmetric SQUID. The position of 
minima, however, slightly deviates from <I> = 0. We have 
assumed that J\ in the left junction is equal to that in 
the right junction. When those are different from each 
other, the Josephson critical current becomes asymmet- 
ric of <I> even in the symmetric SQUID. Thus the small 
asymmetry found in the experiment might be caused by 
the asymmetry of the two tunnel junctions. The degree 
of asymmetry is negligible when the difference of J\ in 
the two junctions is much smaller than themselves. 

In addition to p x + ip yi the sin^) + i sm(p y )£l^ and 
sin^ +p y ) + i sin(p x — Py) symmetries 4 ^ have been pro- 
posed so far. The perturbation expansion! 4 ^ also indi- 
cates a possibility of the chiral p wave symmetry with 
more complicated structure than sin(p x ) + ism(p y ) and 
sin^ + p y ) + ism(p x — p y ). The characteristic behav- 
ior of the Josephson current in the two types of SQUID 
discussed in Sec. II and III are valid not only for the 
isotropic p x + ip y symmetry but also for anisotropic 
sin^) + i sin(p y ) and sin(^ + p y ) + i sin(p x - p y ) sym- 
metries. The latter two gap function have the four- 
hold symmetry in the momentum space consistently 
with a thermal conductivity experiment 46 . In particu- 
lar, sin^ + p y ) + ism(p x — p y ) is close to the chiral / 
wave symmetry because sin^ + p y ) and sin^ — p y ) 
change their sign six times on the Fermi surface. To dis- 
tinguish one symmetry from another, much more detail 
analysis would be required. At present, we propose the 
angle resolved tunneling spectra for this analysis 4 ^. 



VI. CONCLUSION 

We have studied the critical Josephson current in two 
types of SQUID consisting of s wave superconductor and 
Sr2Ru04. In the p x + ip y symmetry, the critical Joseph- 
son current in the corner SQUID becomes the asymmet- 
ric function of $ because the current-phase relation in 
the two junctions relate to each other by tt/2 phase shift 
in the gap function. We also show that the asymmetry 
remains even in the presence of the chiral domain struc- 
tures in Sr2Ru04- Our results well explain the recent 
experimental findings 25 . 



V. DISCUSSION 

We have shown that the asymmetry of the critical cur- 
rent with respect to $ in the corner SQUID is the evi- 
dence of the p x + ip y symmetry when the period of oscil- 
lations is <I>o- A recent experiment shows a large asym- 
metry of the critical current in the corner SQUIE^S 5 -. The 
experiment also shows a minima of the critical current 
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